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CNGS neutrino beam systematics for ^13 
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Energy spectra, intensity and composition of the CERN to Gran Sasso CNGS neutrino beam for — > i/r and 
Vfi — > i^e oscillation searches are presented. The associated beam systematics, which is the major ingredient for 
the 1/^ Ue search sensitivity, are obtained from the study of the previous CERN WANE. 



1. Introduction 

Over the next five years the present generation 
of oscillation experiments at accelerators with 
long-baseline beams, K2K at KEK 1 , MI- 
NOS at NUMI beam of FNAL 2^ and ICARUS 
and OPERA at CNGS CERN to Gran Sasso 
neutrino beams ^ is expected to confirm the 
— > Vr transitions observed in the atmospheric 
1/ and measure sin^ 2^23 and jAmlgj within 10 % 
of accuracy if IATO23I > ^V^. K2K and 

MINOS are looking for neutrino disappearance, 
by measuring the survival probability as a 
function of neutrino energy while ICARUS and 
OPERA will search for evidence of Vr interactions 
in conventional beams. 

Even not explicitly optimized for t'^ 
searches these experiments can explore sin^ 26*13 
oscillation parameter beyond the CHOOZ limit 
0. In particular ICARUS and OPERA can 
reach a combined sensitivity on sin^ 2^13 (CP vi- 
olation and matter effects not accounted for) a 
factor ~ 5 better than CHOOZ in the allowed 
parameter region of atmospheric neutrino oscilla- 
tions for five years exposure to the CNGS beam 
at nominal intensity [Sj. The — ^ Ve oscillations 
will be searched for as Ve charge current events ex- 
cess over the contamination of the beam: the 
key issue will be the knowledge of the neutrino 
beam composition and spectrum. A strong im- 
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pact on sin^ 2^13 sensitivity is expected from the 
systematics of the v^/v^ ratio, and in particular 
from its normalization error. 

The CNGS beam-line design was accomplished 
on the basis of the previous WANE beam ex- 
perience at CERN SPS [T] for CHORUS and NO- 
MAD experiments [S] which allowed for a power- 
ful study of the neutrino beam providing a strong 
benchmark for conventional neutrino beams and 
in particular for the CNGS. In this sense the 
CNGS beam systematics on Vf^jvy, ratio can be 
derived and predicted from the previous measure- 
ments and studies performed with the WANF. 

2. The WANF: a case of study 

In the WANF v^^ facihty 450 GeV/c protons 
were extracted by the SPS and sent into a seg- 
mented target composed by 11 Be rods, 10 cm of 
length and 3 mm diameter each. The produced 
positive (negative) mesons, essentially tt and K, 
were focused (defocused) by two pulsed magnetic 
lenses (horn and reflector) into a 290 m vacuum 
decay tunnel were neutrinos were produced. The 
resulting beam in NOMAD at 840 m from the 
Be target, was characterized by an energy ~ 24.3 
GeV and a contamination of ~ 6.8%, ~ 1% and 
0.3% of F^, and Ve respectively. Four main 
processes as sources of neutrinos were recognized: 

• 450 GeV/c proton interactions in the Be 
target which produce tt^, , K'^,... origi- 
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nating 96% of the i/^ flux in NOMAD; 

• primary protons non interacting or missing 
the target which produce in the horn and 
reflector walls, windows and beam-dump, 
mesons only weakly or not focused by the 
beam-optics 15% of Fj, + T'e); 

• "prompt neutrinos" generated in the decay 
of charmed particles and kaons in the target 
and dump yielding ^ 6% of V^; 



particles reinteractions along the beam-line 
affecting 

defocused component -I- F 



10% of i^f, and i^e but ~ 40% of 
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Figure 2. Neutrino energy spectra (left) for the 
data (points with statistical error bars) and the 
Monte Carlo (histogram) for CC, 17^ CC and 
Ve CC interactions and their corresponding ratios 
(right) in NOMAD 
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relative to the 7r+ is essential to calculate the z^e 
initial content in the beam. 

Direct measurements of tt and K production in Be 
by 400 GeV/c protons were performed by NA20 
Collaboration ^ and by the SPY Collaboration 
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Figure 1. FLUKA to SPY, NA20 meson produc- 
tion corrections as a function of the momentum p. 
The curves are the results of fitting the weights at 
fixed p (showed with their uncertainty bars) with 
combination of polynomial functions 

In order to predict energy spectra, intensity 
and composition of the v beam a sound knowledge 
of meson production in the target is required. In- 
deed the accurate description of the and 
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Figure 3. Energy spectrum for the measured 
h'e CC interactions (points with statistical error 
bars) and the Monte Carlo (histogram) in NO- 
MAD [I3|. 




Figure 4. Schematic layout of the future CNGS neutrino beam hne. 
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[TU] at 450 GeV/c. Accuracy of Monte Carlo gen- 
erators of hadronic interactions for the meson pro- 
duction could limit the sensitivity to neutrino os- 
cillation searches. The best agreement between 
predictions and data was found with FLUKA 
standalone code which reproduces the mea- 
sured TT and K yields at the level of ^ 20% in 
the momentum region 30 100 GeV/c which is 
expected to contribute most to v flux jj^l- Fur- 
thermore an accurate description of the primary 
proton beam spot, focusing system as well as of 
the materials inserted in the beam line from the 
target to the dump (hadronic reinteraction pro- 
cesses) is also mandatory especially for the neu- 
trino beam minor components. 

A complete analysis of this v beam was per- 
formed by a beam line simulation from the Be 
target up to the NOMAD detector including the 
FLUKA generator with further corrections to the 
meson production based on the residual differ- 
ences between the predicted and measured meson 
yield in Be. These reweighting functions were 
calculated as a function of the meson momen- 
tum by integrating the particle production over 
10 mrad of angular acceptance and including also 
the transport efficiency as a function of produc- 
tion angle along the beam-line. The Kq reweight- 
ing was obtained from the corresponding for 
using a "quark-counting model" with 15 % of un- 
certainty [T^ . 

The agreement with the measured v interactions 
in NOMAD was at the few percent level (fig. 0) 
with a systematic normalization error of ~ 7% 
on I'^CC/p.o.t. essentially determined by the ac- 
curacy on the TT"*" and production in the Be 
target (3.4%) and on the proton beam spot posi- 
tion on the Be target. The corresponding energy 
dependent error ranged between 2 and 5% in the 
2 < < 100 GeV neutrino energy. Due to cor- 
relations between the sources of the v^j^ and v^. 
fluxes, the normalization uncertainties on Vel^ii 
ratio was smaller the the uncertainties on the in- 
dividual and fluxes: 4.2% ^21- The cor- 
responding energy-dependent uncertainty ranged 
from 4 to 6 %. 



3. The CNGS neutrino beam 

In the CNGS neutrino beam facility (Fig. ^ 
the primary protons at 400 GeV/c of momentum 
with nominal intensity of 4.5 • 10^^ pot/year (pro- 
ton beam shared operations), will be extracted 
from the CERN SPS and sent to a carbon target, 
13 graphite rods of 10 cm length, 5 and 4 mm of 
diameter, the first eight spaced by 9 cm. Simi- 
larly to previous WANF neutrino beam a mag- 
netic horn and a reflector will allow the focusing 
(defocusing) of positive (negative)charged secon- 
daries into a ~ 1 km long decay tunnel where an 
intense z/^ neutrino beam is produced. A large 
graphite and iron dump will absorb the residual 
hadrons at the end on the beam-line. A partic- 
ular attention was devoted to have the full con- 
tainment of the proton beam spot in the target 
and to reduce as well as possible the quantity of 
material in the beam line in order to maximize 
the neutrino flux in the Gran Sasso site. 

A full detailed simulation of the CNGS beam 
line was performed within FLUKA framework. 
At the nominal proton beam intensity 4.5 • 
10^^ pot/year (proton beam shared operations) 
roughly 2900 z^^CC/kt/year are expected at the 
Gran Sasso site. The increase the proton beam 
intensity has been studied, resulting in a possible 
upgrade bigger than 40% [Tlj . 

The muon neutrino flux will be characterized 
by an average energy of 17.4 GeV and ^ 0.6% Ve 
to contamination for Ei, < 40 GeV (Fig. [SJ. 
The 77^ and 17^ component are below 2% and 0.2% 
respectively. The Vt intrinsic level in the beam 
will be below 10^^ allowing for clean Vr 
appearance experiments with both ICARUS and 
OPERA detectors. Due to the 732 Km of base- 
line the contribution to neutrino beam from the 
is reduced by a factor 1.5^2 with respect 
to the WANF, while those from high energy pro- 
ton interactions downstream of the graphite tar- 
get and from prompt charmed particles and kaons 
decay in the target and dump are negligible. The 
i^e component will be mainly produced in the ^'^ 
decay instead from K'^ and K'^ as in the WANF. 
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4. CNGS beam systematics 

The previous scheme used for the WANF neu- 
trino beam calculation based on FLUKA and 
SPY, NA20 hadroproduction data in Be is well 
suitable for the CNGS neutrino neutrino beam. 
The Be-C scaling due to the different target ma- 
terial is expected to be almost independent on the 
secondary meson type, to depend only weakly on 
the transverse momentum of the meson, so con- 
tributing less than 2% to the uncertainty on neu- 
trino flux and less than 0.5% to Ve jvy^ ratio. 

The systematics of neutrino flux at Gran Sasso 
can be evaluated from the WANF by properly 
rescaling the contribution of each meson to neu- 
trino flux according to the different production 
processes. A total normalization error of ~ 4% 
can be predicted for both v^, and Vf. where ^ 3.8% 
of uncertainty is the secondary production at the 
target and ~ 1.1% error is due to particle reinter- 
actions along the beam-line. Owing to the tar- 
get conflguration and the 732 km base-line, the 
position of proton beam on the target section, 
the beam-line optics such as horn and reflector 
currents and the amount of material included in 
the simulation are not critical but contribute only 
0.8% to normalization error of v^^ and flux. 
Due to correlations between the sources of the 
and fluxes, where both neutrinos are essentially 
produced from tt"*" of different momenta, the v^^ 



by pion decay and v,, from the subsequent muon 
decay, a large cancellation occurs in the v^jv^ er- 
ror evaluation. Therefore the normalization un- 
certainty on Ui-lvfi is expected to be only ^ 3.1%. 
The corresponding bin-to-bin energy dependent 
error is estimated to be 2.5 -j- 4% in the relevant 
part of the spectrum for u^, and fe/fc'/n ratio. 
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Figure 6. Energy dependent uncertainty for i/^ 
and Ve neutrino flux and for their ratio Ue/^n- 
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5. Conclusions 

The CERN to Gran Sasso CNGS neutrino 
beam project for i/^ —> oscillation search with 
ICARUS and OPERA detectors largely benefits 
of the WANE experience with CHORUS and NO- 
MAD experiments. 

The neutrino flux can be predicted at Gran Sasso 
within a systematic uncertainty ~ 5%, the Ve/vij. 
ratio will be known with ~ 3.1% normalization 
error and 3 ~ 4% energy dependent error, better 
than quoted by ICARUS and OPERA, opening 
the possibility to search for — > oscillation 
exploring sin'^ 26'i3 beyond the CHOOZ limit. 
However the sensitivity on both — > Vt and 
— > oscillation channels are completely dom- 
inated by the statistics. In 5 years of CNGS stan- 
dard operations, 4.5 • 10^^ pot /year, 12 Vr CC 
events can be recognized in ICARUS and simi- 
larly in OPERA if Am^g ^ 2.6 • 10"^ gy^. In 
the same time an excess of only 21 — > i^e CC 
events over 120 I'eCC from beam contamination 
and Vr is expected in ICARUS (2.35 kton 

fiducial mass) for E GeV if sin^ IQ^-i ~ 0.04 
(013 ~ 6"). The foreseen increase of proton beam 
intensity is mandatory. 
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